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The historical development is described of the idea that sufficiently strong radio waves 
can modify the properties of the ionospheric regions which they traverse. 

The modified primary properties are the electron temperature and collision frequency, 
the rates of electron attachment and recombination, and the excitation and ionization of 
molecules. 

The resulting changes in the collision frequency and electron density modify the region's 
absorption coefficient and refractive index for any wave, including the modifying wave 
itself. These lead to the phenomena of wave interaction and self -distortion, respectively. 

The associated mathematical theory includes the important effects on gyrowaves 
due to the presence of the earth's magnetic field. The principal predictions of this theory 
and their subsequent experimental verification are discussed. 

Applications of this theory to the study of the ionosphere and laboratory plasmas are 
indicated. These include control of the ionosphere and the generation of artificial airglows 
(or aurorae) by means of powerful gyrowaves. 

The possible modification of the ionosphere by low-frequency extraterrestrial radio 
waves is also considered; such waves include geomagnetic disturbances. They may arise 
from hydromagnetic disturbances in the solar wind and from fluctuating screening by plasma 
clouds of postulated electric charges on the sun. 



1. Introduction 

In a letter to Nature, B. D. H. Tellegen [1933] 
reported that the transmissions of the Swiss radio 
station at Beromunster (650 kc/s) received in 
Holland appeared to be modulated in the ionosphere 
by the radiation from the newly erected powerful 
station at Luxembourg (200 kW, 252.1 kc/s). 
This and an earlier report by Butt [1933] are ex- 
amples of the phenomenon of ionospheric wave 
interaction 2 which was subsequently a source of 
annoyance to European broadcast listeners. The 
observations reported by B. van der Pol [1935] 
and others all indicated that a powerful station of 
sufficiently low frequency can impress its modu- 
lation on the wave received from another station 
(the wanted wave) when the latter traverses iono- 
spheric regions whose projections on the ground 
lie within about 100 to 200 km of the disturbing 
station. The two waves are now commonly called 
the disturbing wave and the wanted wave. The 
observed depth of impressed modulation could 
amount to about 10 percent. 

It is worthy of note by applied physicists that the 
accepted theory of this phenomenon arose ex- 
clusively from a "thought experiment" which 
occurred to the present author while officially 
examining a doctoral candidate's exposition of 
the magneto-ionic theory of ionospheric waves. 
This theory was first mentioned in a letter to Nature 
[Bailey and Martyn, 1934] and more precisely 



developed in subsequent publications [Bailey, 1937a, 
b, c, d; 1938a, b]. A brief account of this develop- 
ment is given in a Supplement to Nuovo Cimento 
[Bailey, 1956]. 

Since most of the experimental work on wave 
interaction has been guided by these mathematical 
developments, it is here desirable to present the 
theory first. 

The two fundamental concepts on which the theory 
is based are as follows: 

1. Any electric field of intensity E(V/cm) applied to 
free electrons in thermal equilibrium with a gas at 
pressure p (mm Hg) and temperature T (°K) will 
impart energy to the electrons and so raise their 
temperature to a higher value T. This results in an 
increase in the rms value u of the random velocities 
of the electrons and so to a change in v the mean 
frequency of collision of an electron with molecules 
of the gas ; in general this change is an increase of the 
collision frequency. 

2. The absorption coefficient k of a radio wave which 
traverses a particular region of the ionosphere varies 
with the value of the local collision frequency v. 
Except for an extraordinary wave of frequency near 
the local gyro frequency, k increases 3 with v. 

It follows from 1 that a powerful radio wave will 
in general increase the collision frequency in some 
region of the ionosphere, and then from 2 that as a 
result another wave (not an extraordinary gyro- 
wave) traversing this region becomes more absorbed 
by the region. 4 This increased absorption of the 



i Deceased. Send communications and reprint requests to J. Bailey, 4 Ch. de 
la Golette, Meyrin, Geneva, Switzerland. 
2 Also known as the "Luxembourg effect" and the "Tellegen effect." 



» For the exceptional extraordinary gyrowave, k decreases as v increases. 
4 An extraordinary gyrowave will become less absorbed by the region. 



wanted wave will reveal itself (in general) as a 
decrease in the amplitude of the wave at reception. 
Thus any modulation of the amplitude of the disturb- 
ing wave at a frequency n and to a depth M will 
impress a modulation of the same fundamental 
frequency on the wanted wave. 

The same concepts 1 and 2 lead to the conclusion 
that a sufficiently powerful wave will modify the 
absorption of the regions which it traverses and so 
be received in a distorted form. This phenomenon 
has been termed "self -distortion." 

Besides explaining the observations of Tellegen, 
Butt, and others, the theory has predicted the follow- 
ing phenomena: 

(1) The amplitude of the wanted wave is modu- 
lated by a fraction j t which is given by formulas 
(44) and (45) and consists of a fundamental part of 
frequency n, phase lag & and modulation depth Mi 
proportional to M and a harmonic part of frequency 
2n } phase lag cp 2 , and depth M 2 proportional to M 2 . 

(2) The depths M x and M 2 of the two impressed 
modulations decrease in specified manners as the 
modulation frequency ?i of the disturbing wave 
increases. 

(3) The depths M x and M 2 increase nearly pro- 
portionally to the mean power P of the disturbing 
station. 

(4) The phase lags lo x and <p 2 are equal to the 
inverse tangents of the ratios n/Gv and 2n/Gv , 
respectively, where v Q is the value of the electron 
collision frequency in the unperturbed local iono- 
spheric region and G is a particular constant derived 
from laboratory data for air; 6 r ^2X10~ 3 . 

(5) The depths M x and M 2 decrease as the angular 
radio frequency co of the disturbing wave increases 
within certain limits. 

(6) Each of the depths M x and M 2 passes through 
one or two maximum values as the angular radio 
frequency co of the disturbing wave is varied about 
the local angular gyrofrequency ft, i.e., resonance 
occurs. Also a radiated power of 1 or 2 kW would 
then suffice to produce observable interaction with 
an appropriate wanted wave. 

(7) The depths M t and M 2 increase as the region 
of ionospheric reflection of the wanted wave ap- 
proaches to about 100 km from the ionospheric 
region vertically above the disturbing station. 

(8) A pulsed disturbing wave causes the amplitude 
of the wanted wave to rise and decay approximately 
according to a particular law which involves the 
time constant t = l/Gv . 

(9) The observed modulation of a wave becomes 
increasingly distorted, at the expense of the lower 
modulation frequencies, as its power increases, i.e., 
self -distort ion increases with the power. 

(10) "Interaction atmospherics" can be impressed 
on a wanted wave by the low-frequency radiation 
from suitably located electric storms. 

(11) Combination tones can be impressed on a 
wanted wave by two disturbing waves with radio 
frequencies which differ by an audio frequency. 



The predictions (1) to (9) have been confirmed by 
subsequent experimental observations. Also (4) 
has been used by Ratcliffe and Shaw [1948] to 
determine the height at which a particular wave 
interaction occurs and (8) has been used to determine 
the quantity Gv which relates in part to the mean 
energy lost by an electron at a collision with an 
ionospheric molecule. The prediction (11) appears 
to account for certain observations of Gill [1935]. 

Also it may be mentioned that the formula for 
M 1 yielded by the theory [and given below under 
(45)] correctly accounts for the order of magnitude 
of the depth of the modulation impressed on Ber- 
omunster's wave by the Luxembourg station as 
observed by van der Pol and van der Mark [1935]. 

The fundamental concepts 1 and 2 have also led to 
the following predictions: 

(12) A specified vertical beam of extraordinary 
gyrowaves with radiated power of 500 kW can 
generate a visible glow discharge or artificial "aurora" 
in the ^-region at night [Bailey, 1938a, b; 1961]. 

(13) A vertical beam of extraordinary gyrowaves 
radiated with a power of 500 kW from an array of 
40 horizontal dipoles can penetrate into the nocturnal 
E'-region much more easily than a weak beam. 
Also, this beam suffices to increase the electron 
density near the 92 km level of the nocturnal E- 
region by an observable amount and thus to deter- 
mine the frequency of attachment of electrons to 
molecules near this level [Bailey, 1959]. 

These results and the corresponding experimental 
studies by Goldstein et al. [1953a, b], and Anderson 
and Goldstein [1955], of the effects of strong micro- 
waves on gaseous plasma have led to a proposal for 
control of the ionosphere by means of radio waves 
[Bailey and Goldstein, 1958]. 

We shall now consider in detail according to the 
following plan the principal theoretical and experi- 
mental results of investigations on ionospheric wave 
interaction published up to February 1963. 

Section 2 summarizes the relevant concepts, 
methods, and experimental observations of the 
Townsend School, which relate to the collisions of 
electrons w T ith gas molecules or atoms, under the 
influence of electric fields of effective intensity 
E V/cm at gas pressures p mm Hg and gas tempera- 
tures T (usually about 15 °C). 

Section 3 summarizes the relevant concepts and 
experimental observations, of the Apple ton School 
and other ionospheric physicists, which relate to the 
propagation of radio waves in the ionosphere regarded 
as a magneto-ionic gaseous medium. 

Section 4 briefly sets out the mathematical theory 
of wave interaction as a logical consequence of the 
work described in sections 1 to 3. 

Section 5 considers in detail the principal predic- 
tions to which this theory leads. 

Section 6 describes the experimental tests and 
verifications of these predictions. 

Sections 7 and 8 discuss the use of the theory as an 
instrument for the study of the ionosphere and elec- 
trical gas discharges, some of its practical applica- 
tions, and some recent developments of the subject. 
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2. Collisions of Electrons With Gas Mole- 
cules Under the Influence of an Electric 
Field 

In the absence of an electric field, the mean energy 
U of an electron in a gas at uniform temperature T 
is equal to the mean energy U of a molecule 5 of the 
gas. 

Townsend and his collaborators [Townsend, 1947; 
Healey and Reed, 1940] discovered that in the pres- 
ence of an electric field E, whether constant or alter- 
nating [Townsend, 1932], the mean energy of random 
motion U of an electron is greater than U , and that 
for a given gas, the ratio 

k=U/U 

is a function of the ratio E/p alone, where p is the gas 
pressure. 

By measuring, in a constant field E, the divergence 
in the gas of a stream of electrons which passes 
through a slit, they were able to determine the values 
of the Townsend factor k corresponding to different 
values of E/p. Then by deflecting the same stream 
through a known angle by means of a transverse 
magnetic field, they were able to measure the corre- 
sponding values of the drift velocity W of the elec- 
trons in the direction of the electric field. 

From these known values of k and W and by 
means of simple 6 approximate theoretical formulae, 
they then calculated the mean free path L of an 
electron of energy U in the gas at the standard 
pressure of 1 mm Hg and the mean proportion \ of 
its own energy which such an electron loses at a 
collision with a molecule. 

From the value of L corresponding to a given 
value of U, the mass m of an electron, and the rela- 
tion U=%mu 2 , we can then calculate in succession 
the corresponding mean random velocity u and colli- 
sion frequency v at a pressure p by means of the 
formulae 

u=(2U/m)i=(2kU /m)t (1) 



v=u/l=up/L, 



(2) 



since the mean free path I at the pressure p is given 

by 

l=L/ P . 

Also, from the corresponding value of X we can 
calculate rj, the mean energy lost by an electron at a 
collision, by means of the formula 



:\U. 



(3) 7 



5 In this paper the term "molecule" is intended to include "atom." 

e The more complicated formulae proposed later hy various workers have not 

heen found necessary in the present state of experimental knowledge of wave 

interaction. 
7 1\ can also be obtained directly from the drift velocity W by means of the 

formula 7^ = 1. 23 mW2. 



The values of k, W, L, and X corresponding to 
different values of E/p for well-dried air were first 
determined by Townsend and Tizard [1913]. These 
are also set out in the books by Townsend [1947] 
and Healey and Reed [1940]. More recent deter- 
minations of some of these quantities have been 
made by Bailey [1925], Neilsen and Bradbury [1937], 
and Crompton, Huxley, and Sutton [1953]. On ac- 
count of the limitations set by the experimental 
methods used, the lowest value of the Townsend 
factor measured in air was Z;=4; this corresponds to 
a mean electron energy of 4/27 electron volts. Since 
the variations of L and rj with U for k<^4 are not 
known, the application of the known experimental 
data to the lower ionosphere makes it necessary to 
extrapolate these data in some simple way to the 
range in which 1<&<4. For this purpose the sim- 
plest two hypotheses to adopt are: 

(1 ) L is independent of U and 

(2) v =G(U-Uo), (4) 

where G is a constant determined by the data ob- 
tained with the lowest values of the field E. For 
air, the value 6r^2X10~~ 3 has proved the most 
successful. 

From the first hypothesis and (1) and (2), it follows 
that the collision frequency v in air always increase* 
with the intensity of the field E in the range of 
energies considered. 

Following Townsend's work on the uniform column 
in high frequency electrical discharges through gases 
(''Electrical Discharges, " Llewellyn Jones [1953]) 
the formulae (1) to (4) may be applied to the effects 
caused by a strong radio wave in the ionosphere. 
But the dependence of these effects on the field 
intensity E, the radio frequency oo/2ir, and the 
terrestrial magnetic field H needs further considera- 
tion, such as that given below in section 4 and the 
appendix. This involves the determination of w, 
the mean work done on an electron over a free path 
or, alternatively, the mean power a=Nvw expended 
on the N electrons in unit volume. 

3. Wave Propagation in a Gaseous 
Magneto-Ionic Medium 

In Appleton's well-known magneto-ionic theory 
of wave propagation, the medium is a gas which 
contains N electrons/cm 3 and is pervaded by a 
uniform, constant magnetic field H; also the mean 
collision frequency v of an electron in the gas is 
taken as independent of the amplitude of the wave. 

The complex refractive index M—\x—ix and the 
polarization number K, which define a mode of prop- 
agation of a wave of angular frequency co in this 
medium, are determined by means of Maxwell's 
field equations combined with the equation of trans- 
port of mean momentum, which governs the momen- 
tum Nm\ of the electrons in unit volume. When 
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N is uniform and constant, the latter reduces to the 
equation 

m -Y-+v7n\=eE+- vXH, (5) 

at c 

where e is the charge on an electron in e.s. units and 
the magnetic vector of the wave field is neglected 
by comparison with H. 

From the dispersion equation yielded by (5) and 
the field equations, it follows that there are two 
wave modes and that for each the corresponding- 
ray becomes attenuated as it progresses, through a 
factor a s defined by 



= exp 



(-/» 



(6) 



where s is the length of the ray path measured from 
its beginning and k is related to x as follows : 



--0)X/C. 



(7) 



The value of the coefficient of absorption k 
depends on v, N, and the angular gyrofrequency 12. 
The theory shows that when 



2 » 2 



(8) 



k is in general approximately proportional to the 
product vN. For the exceptional case, in which 
the angular wave frequency co is nearly equal to il 
and the mode is extraordinary, k is given by (64) 
and so is nearly proportional to iW/{V+(co— ft) 2 } . 
It follows that with increasing v the attenuation in 
general increases, but decreases for an extraordinary 
wave mode with a frequency near the gyrofrequency. 
The latter may be conveniently termed an extraor- 
dinary gyrowave. 

In a slowly varying nonuniform medium a ray 
theory of propagation may be used. Then along a 
ray path of length s, the field amplitude is given by 



where 



E=E exp (—a) 
a= ads. 



Since the mean power flux density is P where 



4. Mathematical Theory of Ionospheric 
Wave Interaction 

From consideration of the results mentioned in 
section 2, we see that in (5) the quantity v is not 
really a constant, but must now be regarded as a 
function of the wave field intensity. 

It follows that the classical magneto-ionic for- 
mulae for M and R must now be regarded as first 
approximations which are accurate only when the 
waves concerned are weak enough to make v closely 
equal to the normal collision frequency v which 
exists in the unperturbed medium. 

The fact that at least k is a function of the field 
intensity E is sometimes referred to as the "non- 
linearity" of the medium. 8 

In order to obtain higher approximations to M 
and R it is in general necessary to apply lengthy 
perturbation procedures. So in order to keep the 
present mathematical theory sufficiently simple for 
practical purposes we shall relax its rigor to the 
limits allowed by the uncertainties in the physical 
data involved. 

4.1. Work w Done on an Electron by a Plane- 
Polarized Wave Over a Mean Free Path 

It has been shown [Bailey, 1937b] that 



w= 



2m 



E 2 e [A+A cos 2co£-£ sin 2a>t] (11) 



where E e is the effective value E/-yj2 of the local 
electric field intensity E, 



A 

B 



2c 2 



*/ 2 +(co-£2) 2 ' *> 2 +(co+S2) 2 ' v 2 +u 2 



/ co— 12 

V+(co- 



co+12 



SI) 2 ' *> 2 + 0+ft) 2 7 v » 2 +co 2 v 



) 



2co c 2 



(12) 



co is the angular frequency of the wave, 
ti=He/mc (angular gyrofrequency), 
s=sin yp, c=cos \p 

and \f/ is the angle__between E and H. 

The mean work w done during a complete cycle 
27r/co and over a free path is therefore given by 



P=c/J5?/4tt, 



w=^ E\A. 

2m e 



(13) 



where E T is the rms value of the component of E 
transverse to the direction of propagation, it follows 
that 

P/n=(Polno) exp (-2a) (9) 

and so 



b(P/M)/a*=-2i(p/ M ). 



(10) 



In the following sections we shall assume that 
E e varies appreciably only during a period of time 
which is large compared with the period 27r/co of 
the wave. 



s The experiments of King, described in section 6, show that notable iono- 
spheric self-distortion can occur when the radiated power of the transmitter is 
at least 200 kW. 
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4.2. Differential Equation for the Mean Kinetic 
Energy U of an Electron 

The equation of transport of energy in the electron 
gas is 

d(NU) . AT AT 

so when the wave is not strong enough to change N 
appreciably we have 



dU 
dt 



-\-VTl = VW. 



(14) 



When the energy U does not exceed 4Z7 we have 
1<&<4 and so we may adopt the formula (4) as a 
working hypothesis. Then (14) leads to the equation 



dU 



+ (p/p )(U-Uo) = tovW 



where 



U=VG»o\ 

T=t/t J" 



(15) 
(16) 



t is a time constant and r is the time measured in 
units equal to this time constant. 
Also, by (2) and (1) 



u/v =(UIUo)KU/L). 



(17) 



The dependence of L on U for [7<^4[7 is not yet 
known with certainty. But the experiments on air 
of Townsend and Tizard and of Crompton, Huxley, 
and Sutton indicate that no great error is likely to 
result from our adopting here the original hypothesis 
that L is independent of U, i.e., 



v/vo=(P/Uo)*. 



(17.1) 



By substitution for v from (17) or (17.1) in (12) 
and from (11) we can express vw as a known function 
of U and E 2 . Then (15) becomes a differential 
equation in the unknown function U with E e as a 
given slowly varying function of t, of the form 



E e =E el M(t)=E el M(t r) 



(18) 



where E eX is a constant field. M(t) may be termed 
the modulation function or modulation envelope. 

The resultant differential equation requires some 
method of approximation for its solution. From 
such a solution U and (17) or (17.1) we can then 
obtain an approximate expression for ^asa function 
of t, F e i, co, 12, \[/, G, and U or v . 

These expressions for U and v contain terms 
which vary slowly with t and terms which are 
periodic in t with 2co as their fundamental frequency. 
Retention of the latter terms in the discussion of 
wave interaction and self -distortion would lead to 
similar periodic terms in these phenomena. But at 
radio frequencies co such terms are usually negligible 
by comparison with the slowly varying terms, so 
henceforth they will be neglected. 



Accordingly we will now consider only the varia- 
tions of U and v which take place in a time dt large 
compared with the wave period 27r/co. This requires 
the symbol w in (14) and (15) to be replaced always 
by its mean value w which is given by (13) or by (48). 

4.3. Differential Equation for the Mean Collision 
Frequency v 

It is now convenient to base the discussion on the 
equation for v. Thus, on eliminating U from (15) 
by means of (17.1), replacing w by w, and using (13) 
and (18), we obtain 



^+l(a 2 -l)=w/2GU =RM 2 (t r) 



where 



v/vq, 



R- 



~4GmU ( 



E 2 el A. 



(19) 

(20) 
(21) 



Here a is the collision frequency expressed in units 
equal to v . 

On account of its factor A, the quantity R in (19) 
is a function of a 2 such that this equation has no 
obvious solution in finite terms. 

However, in order to apply the theory to most of 
the observations of wave interaction which have 
been made hitherto, it is sufficient to consider only 
those approximate solutions which correspond to 
values of w/2GU or RM 2 (t T) such that a varies 
about some constant value t^ by a small fraction of 
<ri, i.e., 



where 
Thus 



<T = <Ti + 8, 

v = vi + v 8, 



|«|«<Tl 

M|<<>i. 



(22.1) 



(22.2) 



These approximate solutions are then found from 
solutions of the equation: 



where 



f T + a l 5=R l M\Ur)+\{\-a\) 



Ri ~46mU E2elAl 



and 



A r - 



2c 2 



'vi+(a>-ny- ■ vi+iu+a) 2 ' vi+c 



(23) 



(21.1) 



(24) 



4.4. Collision Frequency v s in the Steady State 
Under the Influence of a Continuous Wave 

Here we have M 2 (t T) = 1. Then the steady value 
ai = Vs /v is obtained by setting da/dt=0 in (19) and 
solving the resultant quartic in o 2 or by setting 5=0 
in (23) and solving the same quartic in a 2 .. 
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When E el is not too large, 9 we then deduce the for the steady state corresponding to the mean 
following approximate expression for the steady effective force 
collision frequency: E el (l + ^M 2 )* 



where 



v s = v +R v , 



En 



4:GmU 



E 2 el A 



(25) then the constant part of the R.H.S. of (23) vanishes, 

Bi(i+W 2 )+Ki-^)=o, 



(26) 



and so (23) reduces to 
d8 



and A is the same as Ax in (24) with v\ replaced by v Q . 

From (25), (26), and (24) it follows that, except 
when the wave's field vector E is only slightly 
inclined to the constant magnetic field vector H, 
v s passes through a maximum value as the wave 
frequency co is varied around the gyro frequency £2. _, , , ,^ 

The corresponding mean electron energy U s behaves a > ^ \^5), 
similarly. 

These resonance effects and their consequences 
are of great importance, as will be shown below. 



-j^+(T 1 d=2R 1 M sin {nUr) — \EM 2 cos (2^ r), 
where 

E 2 el (i+m 2 )Av Q . 



(29.2) 



(30) 



Vs = V0-± 



4GmU 



(31) 



4.5. Rise and Decay of v Caused by a Square Pulse 
of Radio Waves 

When v=v at £=0 and M 2 (t r) = l for 0<t<T 
then (23) has oi=l and its right-hand side constant. 

The corresponding solution of (23) which describes 
the rise of v, is 



The general solution of (30) combined with (22.2) 
leads to the following formula for v\ 



, BMsm. (nt—pi) jBM 2 cos (2nt—<p 2 ) 



(l+* 2 tf)* 



(l+4ra 2 tf)* 



where 



v=v s —{v s 



-v^e-^K 



(27) 



B 



e 2 E 2 el A s t s 



Similarly we obtain the following solution which 
describes the decay of v from the steady value v s : 



where 



v = v +(v s — v Q )e tft s 



(28) 



m 2 l 2 

t s =l/Gv s 

tan (pi = nt s , tan <p 2 =2nt S)> 



(32) 
(33) 



A s is the same as Ax in (24) with v x replaced by v s 
and C is a constant determined by the condition 
v=vq at t =0. 

For the E'-region t is of the order of 1 msec or less. 
Hence the rise and decay require about 1 msec to be 4 - 7 - Modulation of a Wanted Wave Through the 



practically completed. 

Accurate expressions for the rise and decay can 
also be obtained, since (19) is now exactly integrable 
[Bailey, 1937b, eq (34)]. 

4.6. Modulation of v Caused by a Wave Which Is 
Sinoidally Modulated in Amplitude 



Variation of the Local Collision Frequency 

In ^ passing through the ionosphere, the field 
amplitude F of the wanted wave mode is attenuated 
according to the formula 



where 



For such a wave we have 

M(t) = 1 +M sin nt 
0<M<1, n«<a 



F=F exp (-a) 

a= K(v)ds, 
Jo 



(34) 



where 
Then 



M 2 (£ r) = l + iM 2 +2Msin (nt r) 

— JM 2 cos (2nt r). 
If we take 



(29) 



(29.1) 



s is the length of the wanted wave's path measured 
from its entry into the ionosphere, F is the field 
amplitude at entry, and k(v) is the absorption 
coefficient corresponding to the value v of the 
collision frequency at the point s in the path. 

Under those conditions in which the collision 
frequency varies about some constant value v u i.e., 



we have 



9 More definitely when i?o«^- 



a=ai-\-Aa 



(35) 
(36) 
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where 



and 



<xi— I n{v\)ds 

Jo 

Aa= I {k(v) — k(vi) }ds 
Jo 



(37) 



We now apply (41) to determine the modulation 
fractions j Y and f 2 of the wanted wave respectively 
during and after the disturbing square pulse con- 
sidered in subsection 4.5. 

Since v=v Q at t=0 and M 2 (t r) = l for 0<t<T, 
we then have <r 1 = \ and so (41) yields the formula 



Hence, by (34) and (35) 

F=F X exp (— Aa) 
F^Foexp (-a,). 



/!=&(!-«-"'•). 



(43.1) 



where 



Similarly, if v=v s at t=0 and M 2 (t T)=0 for 
ty>T, we have ai=v s /i> and so (41) yields 



f2 — Si 



(43.2) 



So at the point s the amplitude of the wanted wave 
is modulated according to the modulation function 
M t (t) where 



If we set 



Mi(t)=exp (— Aa). 
Af,(0=l+/i(0 (38) 

the impressed modulation fraction /<(£) is given by 

/z(0=--exp(-Aa)-l. 



Next we apply (41) to determine the modulation 
fraction j i impressed on the wanted wave by the 
sinoidally modulated disturbing wave considered in 
subsection 4.6. 

As in 4.6, we now have: M 2 (t r) given by (29.1), 



i.e. 



<ri = <r s =Vs/vo, 



and (29.2) which is satisfied by the root <ti = <t s 
From (42)/ (29.2), and (29.1) we obtain 



When observation shows that[/,(*) |«1, we will S l M 2 (t r) + S 2 =S 1 {2M sin (nt r)-W 2 cos (2nt Q r)} } 

necessarily have the approximation 

and so (41) now becomes 

^+aJ=2S l Msin (^ r)-^iM 2 cos (2nt Q r). (41.1) 

The impressed modulation fraction is given by the 
general solution of this equation; namely, 



/,(«) = -A*= - P {K(p)-K(v 1 )}ds. (39) 

On substituting for v from (35) and expanding 
k{v x -\-Av) ) we obtain 



/«(*)= -J *k'W Ai 



ds 



provided that 



\Avk"M/k'M\«2- 



When RM 2 (t) is sufficiently small for (22.2) to 
hold true, we have, by (35), 



(40) j t =M x sin (nt—<pi)—Mi cos (pnt-^ + Ce-'"* 
where tp Xi p 2 , t s are defined under (33), 

M l =BMKl+n'tl)*, 

M 2 =iB s M 2 /(l+4nH 2 8 )i 

B s =—2<Ti 1 Si 



(44) 

(45) 
(46) 



where 8 satisfies the equation (23). So, on multiply- and C is a constant such that /,=0 when t=0. 
ing (23) throughout by -k'Mv^s, integrating Since Vl = v, we find from (21.1) and (46) that m a 

with respect to s, and using (40) we obtain the re S lon in whlch "o and E { are constant 
following equation for /=/ ( (i), for a situation in 



which <ti is nearly constant throughout the region 
concerned : 



B. 



e 2 E 2 el A s 



I w n'{v s )ds w 



(46.1) 



dj 



+ c l j=S l M 2 {t r) + S 2 



(41) 



where 



Si = — I Rik' '(vi)v ds 
S 2 =\£(al-l)K>(v l ) 



Vtf 



ds 



where A s is the same as A { in (24) with v s replacing 
vi and s w is the length of the path of the wanted wave. 

It will be noted that all the quantities in (44) which 
relate to the wanted wave are contained in the factor 
B s which occurs in the quantities M x and M 2 . 

A simple calculation shows that when, for the 
( ) wanted wave, k{v) is approximately proportional to 
v the formula (44) confirms a previously given 
formula for M x [Bailey, 1937a, part 1, eq (42)]. 
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4.8. Modulation Impressed by a Wave Mode 

When it is necessary to take into account the 
attenuation of each of the wave modes into which a 
plane wave divides, the theory given above must be 
modified by using in the place of the formula (13) 
for w a formula for either mode which is obtained as 
follows [Bailey, 1937a, part 1]: 

The mean power E supplied by the wave to the 
N electrons in unit volume is given by 






(47) 



where P is the mean power flux density and s d is 
the length of the path of the disturbing wave from 
its place of entry. 

Since this must be equal to N times the mean 
power vw supplied to one electron, it follows that 



Nvw='E. 



(48) 



Also, when fx varies slowly with s d we have from 
(10) and (9) 



E=2jfcP=2jfcPo exp (-2a d ) 



where the subscript d to a symbol relates it to the 
disturbing wave. 

Since the amplitude of the wave mode is modulated 
by the factor M(t) as in (18), therefore 

P Q =P Q1 M 2 (t) 

where P i is the mean flux density of this mode, over 
one period, at entry into the ionosphere. 
Hence, by (49) 



2=3^(0 



(50) 



where 



Ei=2/c rf Poi exp (—2 KadsA (51) 

and K d = K d (v) when expressed as a function of v. 

On using the value of w given by (48) and also 
(50) the equation (19) assumes the form 



where 



f r +l {o*-\)=RM\Ut) 



<t=v/v , T=t/t =Ov t, 
R= l S i /26U N Vo cT. 



(52) 

(53) 
(54) 



As in subsection 4.3 we shall be content here to 
consider only those situations in which R M 2 (t T) is 
small enough for the condition (22.1) to hold true. 



Then (52) may be replaced (approximately) by the 
equation 



dd 



where 



j- T +a 1 8=R 1 M 2 (t r)+ f(l-*?) 



B l =S l (v 1 )f2GUoNv 1 



(55) 



(56) 



and Ei(^i) is the value of Si in (51) when K d =n d (vi) 
and vi=v . 

Since (55) has the same form as (23) we can now 
proceed, as was done with the latter, to derive from 
it the equation for the modulation fraction j=j t {t) 
impressed by the wave mode under consideration 
when the ionospheric region concerned is nearly 
uniform in temperature and pressure and o x is nearly 
constant throughout this region. 

This procedure yields (41) for / and the formulae 
(42) but with Rx now given by (56); i.e., S x is now 
given by the formula 



(49) where 



S ^'2mS^ dS (57) 

q = K , M/Nv l (58) 



and Ei is given approximately by (51) with v = v$. 

We thus see that the formulae (43.1), (43.2), and 
(44) for the impressed modulations j u f 2 , /*, respec- 
tively, are also valid here when S x is given by (57). 

The only factor in S x which relates to the dis- 
turbing wave is Si. For this reason it has been 
termed the index oj interaction. Si also represents 
the mean input of power per unit volume of the iono- 
sphere per period of the disturbing wave. 

The combined effects of the two modes which con- 
stitute a plane wave can be similarly calculated if 
the distribution of the energy flux between the 
modes is known. 

We now consider the common case (A) in which 
the wanted wave is an ordinary mode with 



k(v) = pNv 



(A) 



where p is nearly constant throughout the region 
concerned. In this case we generally have 

co 2 >>. 2 and w 2 >>47rA r 6 2 /m. 

Therefore, by (58), q = p/vi and so q is now nearly 
constant throughout the region. Hence, by (57) 
we now have 



S 1 = 



Sic 

1 Jo 



(57.1) 



2GUoN Vl . 
Since, by (50) and (47), 

M\t)\ H,(fe=- secjSv-cfe, (59) 

Jo Jo Os d 
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where is the angle between the elements of path 
ds d and ds contained between adjacent wave fronts 
of the disturbing wave, 10 we cannot evaluate the 
integral in (57.1) in finite terms except when fi is 
constant (or nearly so) throughout the region of 
interaction; this will occur in the following sub- 
cases, Al and A2. 

Al. Both waves are propagated vertically and 
then |sec/3| = l throughout. 

A2. The disturbing wave is propagated vertically 
and the region of interaction occurs mainly where 
the wanted wave is undeviated. 

In these cases (59) yields 

M 2 (t) £Ztds=sec /3(P -/\) 

where P and P Sd are the values of the disturbing- 
wave's flux density at the beginning and end of a 
length s d of its path. 
Since 

P =P 0l M 2 (t) and P 8 =P 8dl M 2 (t) 

where P i and P , v i are the mean flux densities over 
one period at, respectively, entry into and exit from 
the ionosphere, it follows that 



It follows from (46.1) that B s behaves similarly, 
and that in consequence each of the modulation 
depths M x and M 2 passes through a sharp maximum. 
By (21.1), (42)," and (43.1), the modulation j x 
impressed by a pulsed gyro wave will behave similarly. 

Since the two component modes of the plane 
gyro wave are absorbed at very different rates as they 
propagate through the ionosphere, a more precise 
discussion requires the effects caused by each of these 
modes to be considered separately. We may then 
use the results obtained in subsection 4.8. 

For our present purpose it is sufficient to consider 
only the effects caused by the much more absorbable 
extraordinary component u in the lower E region of 
the nocturnal ionosphere, at moderate latitudes, for 
which we have 



A^<800, G = 10 7 , *><10 6 . 



(63) 



It has been shown [Bailey, 1938a] that the refrac- 
tive index n of the extraordinary mode defined by (61) 
is closely equal to unity and that its coefficient of 
absorption K d is given approximately by 



where 



x 



_ K(l + cos 2 d)Nu 



#=0.0266, 



(64) 



Ei«fe=sec 0(Poi-P, d i), 



and so from (57.1) that, in the subcases Al and A2, 

(60) 



„ _^_ k(v ) sec (T) - 

bl ^~26U N Vl U 01 ~ 7 s ^' 



and 6 is the angle between the directions of propaga- 
tion and of the earth's magnetic field. 

More recently K. A. Smith [1957] has obtained 
the following closer approximation: 

When 



N<500, 12 = 10% 3X10 5 O<1.2X10°, !#!<65 

then 

K(l+cos 2 $)Nv 



Under A2 we have sec /3=sec i w where i w is the 
angle of incidence of the wanted wave. We then 
obtain from (44) and (60) a formula for the depth 
Mi, of the impressed fundamental modulation, which 
is equivalent to the equation (28) given in Huxley's 
Synopsis [1952] and is originally due to Ratcliffe where 
and Shaw [1948]. 

4.9. Gyro Interaction: Interaction Caused by an 
Approximate Gyrowave 

When the disturbing wave is plane-polarized and 

has a frequency within 20 percent of the local gyro- F rom (64) or (65) it is clear that under the specified 

frequency, i.e., conditions K d passes through a notable maximum 

co=i2(l + e) ^ value as co varies about the value Q. On the other 

where ^, (61) hand (51) shows that Hi may then pass through 



p 2 + e 2 ir 2 -2f (A l e-A 2 e 2 ) 

$=3.19X10W 

A=i(2+3sin 2 , 
A 2 =h(l+3 sin 2 6) j 



V- 



(65) 



(66) 



<0.2 



)- 



the formula (24) with vi = v s shows that in a region 
where 

v]«& (62) 

the quantity A s passes through a sharp maximum as 
the wave's frequency co is varied about the gyro- 
frequency ^. 

io For which ds d =ds cos 0. 



either one or two maxima according as the path 
length s d is shorter or longer than some particular 
value. As shown by (57) this behavior will be 
reflected in S x through the integral 



Ziqds, 



11 When required, the weaker effects due to the ordinary component can be 
calculated in a similar way. 
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and therefore, by (43.1) and (44), the modulations/ 
impressed by a pulsed or sinoidally modulated gyro- 
wave will pass through one or two maxima as the 
disturbing wave's frequency varies about the gyro- 
frequency. The curves representing Si, M A , or M 2 
in terms of e, co, or co/27r, may therefore be called 
resonance curves. 

Whether the extraordinary gyrowave is trans- 
mitted from the ground or from a rocket or a satellite 
in (or above) the ionosphere, calculation shows that 
in general the wave is almost completely absorbed 
within a region of the ionosphere about 6 km long. 

Consequently we may expect the interaction to be 
greatest when the ionospheric path of the wanted 
wave lies within this region. 

With the gyrowave transmitted from the ground 
this region at night lies in a horizontal slab with its 
lower face at about 90 km and, for e=0, its upper 
face at about 95 km above the ground. As | e \ 
increases from to 0.2 the upper face rises from 
about 95 km to about 97 km. 

Consequently, when the top of the wanted wave's 
path lies below 95 km, the resonance curve for the 
impressed modulation in terms of e possesses only 
one maximum (dromedarian resonance) . But when 
this top lies above 95 km the resonance curve 
possesses two maxima nearly symmetrically situated 
about the minimum which occurs near e=0 (bactrian 
resonance) , for increase of | e | from initially causes 
the slab to embrace an increasing portion of the 
wanted wave's path and so increases the interaction. 

This deduction is supported by the published 
curves for the index of interaction Si at different 
levels in a specified model of the lower ionosphere 
[Bailey, 1 938a] . For these have one maximum at the 
lower levels and two maxima at the higher levels. 

For a particular model of the nocturnal ionosphere 
and the conditions relating to experiments on gyro 
interaction at Armidale, New South Wales, during 
the period 1950-1954, R. A. Smith [1957] has carried 
out the computations which are needed to evaluate 
qds and the integral in (57). These computations 
led to a theoretical resonance curve for the impressed 
modulation which was in agreement with the 
observed resonance curves [Bailey, Smith, et al., 
1952] within the limits of experimental error. Some 
of these curves are shown in figure 3 below. 

When the disturbing and wanted waves are both 
propagated vertically, the computation of the 
integral in (57) is notably simplified. The experi- 
mental work is also much simplified provided that 
the power is increased sufficiently to compensate for 
the reduction in interaction which would otherwise 
occur. 

If also the wanted wave is an ordinary mode under 
the conditions discussed above, (57) may be replaced 
by (57.1) and so the theoretical computations then 
become simpler still. 

The impressed modulations predetermined by 
means of this theory of gyro interaction have been 
fully confirmed within the limits of experimental 
error by subsequent experiments. 



5. Predictions From the Mathematical 
Theory 

The formulae (44), (45), and (46.1) lead to the 
predictions labeled (1) in section 1. 

The formulae (45) specify the manners in which 
Mx and M 2 decrease with increasing n, as stated 
in the prediction labeled (2). These variations 
with the modulation frequency are similar to the 
variation of the current in a choke coil as the fre- 
quency of the applied p.d. is varied. 

The formula (46.1) shows that when E 2 el is not 
too large, B s is nearly proportional to E 2 el and, 
therefore, also to the mean power P of the disturbing 
station. It then follows from (45) that M x and M 2 
are nearly proportional to P and so we obtain the 
prediction (3). 

The formulae (44) and (33) show that 

<Pi = tan -1 (n/Ov s ), <p 2 = tan -1 (2n/Gv s ) . 

Hence, when the mean power of the disturbing wave 
is such that v SJ the steady value of the perturbed 
collision frequency, differs only slightly from v Qy 
the unperturbed value, we obtain the prediction (4). 

The formulae (45) and (46.1) show that M x 
and M 2 are proportional to A s , where A s , is given 
by (24) with vi replaced by v s . Hence, when v & 
differs only slightly from v , we obtain the predicted 
behavior (5) and also the single maxima predicted 
under the label (6). But the more precise dis- 
cussion of subsection 4.9 shows that when s d the 
path of the disturbing wave exceeds a certain value 
each of the modulations depths passes through two 
maxima. Also, computations by means of the 
formulae given in subsection 4.9 show that when 
€=0 a mean power of 1 or 2 kW radiated from a 
vertical aerial suffices to produce observable inter- 
action. All this is summed up in the prediction (6). 

The rays from an average radio station which 
produce the largest values of E 2 el in the E'-region 
are inclined at 45° to the vertical. Consequently 
the ionospheric regions where a wanted wave is 
most modulated by the disturbing station will be 
about 100 km from the ionospheric region vertically 
above this station. This is the prediction (7). 

The formulae (43.1), (43.2), (16), and (33) show 
that, when the disturbing wave is such that v does 
not differ much from v , the modulation impressed 
by a pulsed disturbing wave rises and falls approxi- 
mately according to the same exponential law with 
the time constant t =l/Gv . This yields the 
prediction (8). 

As the mean power P of a wave increases the 
modulations of its amplitude with the same depth 
M, but different frequencies n will, by (32), cause 
modulations of v of depths which are greatest for 
the lowest frequencies. Consequently the modu- 
lations of the lowest frequencies will be received 
with the smallest amplitudes. This self-distortion 
must increase with P and so we obtain the pre- 
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diction (9). Quantitative discussions of self-dis- 
tortion have been given by Hibberd [1955, 1957]. 

The prediction (10) follows immediately from the 
fact I hat the low-frequency radiation from an electric 
storm reaches the nocturnal E'-region with great 
intensity and so considerably increases the local 
electron collision frequency for periods of the order 
of one millisecond. 

The same general mathematical theory [Bailey, 
1937b] which yielded (11), for a single disturbing 
wave, also yielded a formula [there labeled (35)] for 
the variation of the collision frequency caused by 
two or more disturbing waves. The latter formula 
contained terms representing combination tones and 
so the prediction (11) immediately follows. 

The predicted phenomena (12) and (13) are of 
great interest and possible future importance for the 
study of the ionosphere and ionospheric radio propa- 
gation. Since no experimental tests of the prediction 
( 12) have as yet been carried out, it is sufficient here 
to refer to the original publications in which it is 
discussed in detail [Bailey, 1938, 1959]. The pre- 
dictions (13) are under experimental test by Pro- 
fessor R. A. Smith at the University of New England 
in New South Wales, as indicated in section 8. 

6. Experimental Tests of the Theory and 
Verification of the Predictions 

The principal experiments which served to test 
one or more of these predictions will now be men- 
tioned in their historical order. 

Van der Pol and van der Mark [1935] measured 
the depth of modulation M impressed on .Radio 
Beromunster by Radio Luxembourg at four different 
frequencies of modulation n. Their observations 
are represented by the four dots shown in figure 1. 
The theoretical curve representing M x as a function 
of n, with appropriately selected values of the two 
constants (B S M) and t s in (45), is also shown there. 
The resulting value of 6v = l/t s then yielded a value 
of v Q which was consistent with the values of v pre- 
viously estimated by E. V. Appleton [1932] and 
S. Chapman [1932] from other experiments. This 
constituted the earliest verification of the prediction 
(2) with regard to M x . Also the measured depths of 
modulation were later found to agree in order of 
magnitude with the theoretical value M x given in 
(45), thus verifying a part of the prediction (1). 

Under the guidance of E. V. Appleton [1934] and 
Ralph Stranger [1934], the members of the World 
Radio Research League in Europe made observations 
which roughly verified the prediction (7). 

Baumler and Pfitzer [1935] qualitatively verified 
predictions (2), (5), and (7). 

Extensive experiments which were carried out in 
1937 with the help of the British Broadcasting Cor- 
poration and other European authorities verified 
at least the last part of the prediction (6) [Bailey, 
1937d]. During these experiments, E. W. B. Gill 
of Merton College, Oxford, reported observations 
which appeared to verify prediction (11). 




400 €00 800 

MODULATION FREQUENCY, c/s 



1000 



Figure 1. Modulation frequency, c/s. 
[Bailey, V. A., and D. F. Martyn, Nature 135, 585 (1935).] 

From 1940 onwards M. (Jutolo of Naples and his 
colleagues successfully made extensive observations 
of gyro interaction using disturbing waves with 
powers less than 1 kW. They also succeeded in 
demonstrating the existence of both double-humped 
(bactrian) and single-bumped (dromedarian) reso- 
nance. Thus they verified the whole of the predic- 
tion (G). 

Commencing in 1946 L. G. H. Huxley and J. A. 
Ratcliffe and their respective colleagues carried out 
extensive experiments, both separately and jointly, 
which verified the predictions (1), (2), (3), and (4) 
(excepting the proportionality of M 2 to M 2 ). More- 
over, Huxley, Foster and Newton [1947] observed 
the gyro interaction caused by the 10 kW station 
at Stagshaw; this is consistent with the last part of 
prediction (6). In addition, Ratcliffe and Shaw 
[1948] directly verified the fundamental concept of 
the theory that the amplitude of the wanted wave 
is decreased when the disturbing wave is operating. 
Also Ratcliffe and Shaw successfully applied the 
theoretical formula for the phase lag <fo, given under 
(33), to determine experimentally the height of the 
ionospheric region at which the wave interaction 
took place. Since the heights so determined were 
of the same order as the heights postulated in the 
original theory, we may regard this work as an ex- 
perimental verification of the prediction (4). 

In 1949 experiments on gyro interaction were 
commenced in Eastern Australia by V. A. Bailey 
and his colleagues at the University of Sydney in 
collaboration with R. A. Smith and his colleagues at 
the University of New England at Armidale 
(N.S.W.). These were planned to conform as 
closely as possible to the published requirements of 
the theory [Bailey, 1937c and 1938a]. Accordingly 
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the arrangement of stations shown in figure 2 was 
adopted. In this, the disturbing gyrowave A was 
radiated vertically from a horizontal aerial at Armi- 
dale, which lies at the midpoint of the 740 km line 
that joins Brisbane and Katoomba. The wanted 
wave B was radiated from Brisbane on 590 kc/s 
and received at Katoomba. The local gyrofrequency 
at Armidale was estimated to be 1530 kc/s. 

The gyrowave was radiated 40 times per second 
as a square pulse of length 1 msec and power 36 kW. 
Its frequency/ was kept constant for 2 min at each 
of 8 values lying within ±23 percent of the gyro- 
frequency 1530 kc/s. Each of the cycles 1, 2, and 
3 ran through the middle 6 frequencies and lasted 
17 min. Cycle 4 ran through all the 8 frequencies 
and lasted 23 min. 

In order to provide a standard of comparison for 
the depth of the modulation impressed on the wanted 
wave B, the latter was modulated at its source for 
30 sec in each 2 min period with a standard tone of 
80 c/s to a depth of 5 percent. This depth was chosen 
because computations based on our mathematical 
theory of gyro interaction had indicated in advance 
that the experimental arrangement shown in figure 
2 would impress on B modulation depths of the 
order of 5 percent. This prediction was nicely veri- 
fied by the experimental observations. 

A detailed account of these Australian experiments 
on gyro interaction was published by Bailey, Smith, 
Landecker, Higgs, and Hibberd [1952]. A briefer 
account was published by Bailey in 1956. Examples 
of the experimental resonance curves are shown in 
figure 3. It will thus be seen that these experiments 
not only verified the predictions (6) and (8), but 
also confirmed the theory quantitatively in every 
particular. 

The theory of interaction, including gyro inter- 
action, was also thoroughly verified by the experi- 
ments on microwave interaction in gas-discharge 
tubes of Goldstein, Anderson, and Clarke [1953] and 
of Anderson and Goldstein [1955, 1956]. In some 
of these experiments Goldstein was able to demon- 
strate also the self-distortion of a powerful wave 
[Private communication, 1955], thus verifying the 
prediction (9) for a plasma. 

The prediction (9), that self -distortion in the 
ionosphere increases with the power and does so at 
the expense of the lower modulation frequencies, 
was verified by J. W. King [1957]. 12 In carrying 



i 2 In this publication King inadvertently contradicts the statement on "Iono- 
spheric self-interaction" previously made by Hibberd [1955]; namely: "Its 
existence was first predicted by V. A. Bailey in 1935, . . ." 



out his experiments he was guided by the quanti- 
tative theory of self-distortion developed by F. H. 
Hibberd [1955, 1956]. 

Other tests of the theory arose from the experi- 
mental study of wave interaction by J. A. Fejer 
[1955] in which both the wanted and disturbing 
waves were pulse-modulated. This study demon- 
strated the existence of interaction during the day at 
heights between 70 and 90 km and also that a modi- 
fication of the theory suggested by Huxley was not 
valid. 

The prediction (5) is now verified by the fact that 
during the last 30 years the greatest disturbance by 
interaction backgrounds has been produced by the 
radio stations with the lowest wave frequencies. 

From the whole discussion given above we see that 
experimental observations have verified the 10 
predictions (1) to (9) and (11). The predictions 
(10), (12), and (13) remain to be tested. We may 
therefore conclude that in the form presented here 
the theory of wave interaction is experimentally 
established. 



7. Some Applications to Research on the 
Ionosphere and Gas Discharges and to 
Radio Communication 

The experiments mentioned above by van der Pol 
and van der Mark; Cutolo, Huxley, Ratcliffe, and 
Shaw; Bailey and Smith; and Fejer, combined with 
the theory, have among them served to measure 
the following important quantities relating to the 
ionosphere : 

The electron collision frequencies v, the electron 
densities N, the local magnetic fields H, the constant 
G (relating to the mean loss of energy by an electron 
at a collision with a gas molecule), and the heights at 
which the wave interaction occurs. 

Similar applications to the D-region, have been 
made more recently by Fejer and Vice [1959], 
Bjelland et al. [1959], and Barrington and Thrane 
[1962]. 

The experiments by Goldstein and his associates, 
mentioned above, also demonstrate the value of the 
phenomenon of wave interaction and its theory for 
investigating several fundamental processes in a gas 
discharge. 

As an illustration of its practical value for radio 
communication, it may be mentioned that the author 
has been consulted about it more than once by the 
leading broadcasting organizations in Europe and the 
United States. Also, in the URSI Information 
Bulletin No. 73, on p. 59, it was reported that the 
C.C.I.R. stressed the practical importance of infor- 
mation on nonlinear effects in the ionosphere. 

The present trend towards the use of increasingly 
powerful beams in connection with the reception 
of irregularly scattered waves has also compelled 
attention to the possible occurrence of observable 
wave interaction and other nonlinear phenomena 
at very high frequencies. 
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Figure 3. [Reproduced from Bailey, V. A., R. A. Smith, K. Landecker, A. J. Higgs, 
and F. H. Hibberd, Nature 169, 912 (1952).] 



Lastly, the possibility mentioned by Bailey and 
Goldstein [1958] of controlling the ionosphere in 
various ways by means of powerful radio waves 
has attracted some attention. 

8. Some Recent Developments 

Filling in the gap between the two situations 
corresponding to wave-interaction caused by mod- 
erate powers and ionospheric airglows caused by 
extremely high powers and large aerial arrays, 
an investigation was made [Bailey, 1959] of the 
possible increases of both the collision frequencies 
v and the electron densities N caused by a powerful 
extraordinary circular gyrowave in the nocturnal 



lower ^-region and in the daytime D-region. The 
theory was based on the most reliable available 
experimental data on the behavior of free electrons 
in air with energies up to 1.22 eV. The principal 
process then freshly introduced was the attach- 
ment of electrons to molecules in some collisions. 
It was found that some of the resultant effects, 
of such increases in v and iV, on this wave or on 
other waves are remarkable in magnitude or in 
kind. Thus, when the gyrowave produces near 
the 92 km level at night a reduced equivalent 
electric force (Z/p) equal to 2.5 V/cm per mm Hg 
(measured at 288 °K), it causes notable selfdis- 
tortion and extremely large wave interaction; 
also within 3 msec its penetration into the E'-region 
increases considerably. 
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Other possible observable effects which may be 
caused at night when Zip has a value between 
2.5 and 8 are changes in the visible and infrared 
airglow, in parts of some trails of meteors and 
fast particles of solar or cosmic origin, in the radar 
echoes from such trails and in the local magnetic 
elements. 

The theory also showed that a pulse of gyro- 
waves such as that specified above could be used 
to determine by experiments on the ground the 
rates of attachment of electrons and other electron- 
collision data pertaining to the regions near the 
92 and 75 km levels. 

For these and other reasons, it was therefore 
suggested that experiments be carried out with 
gyrowave pulses beamed from an array of at least 
80 horizontal dipoles connected to a generator of 
at least 500 k W power. 

As a result of this work and of my recommendation 
to the Air Force Cambridge Research Center at 
Bedford (Mass.), the University of New England 
is under contract with this research center to carry 
out a project to " Modify the Ionosphere by means 
of Gyrowaves" using an array of 40 dipoles connected 
to a generator of 500 kW power. The results of 
this project are being reported by R. A. Smith. 

The power density requirements for the excitation 
of an artificial airglow by means of gyrowaves 
[Bailey, 1938b] have been reconsidered [Bailey, 
1961] in the light of later information. One im- 
portant conclusion was that an array of 80 dipoles 
radiating vertically a beam of gyrowaves with a 
mean power of 500 kW would produce an easily 
observable enhancement of the nocturnal airglow 
overhead within a solid angle of 32 square degrees. 
Another conclusion is that from the same array 
pulses of power 3500 kW and length about 0.2 
msec would each produce a flash of airglow that 
is 40 times as bright as the night sky. 

An extensive report on nonlinear phenomena in a 
plasma was published by Ginzburg and Gurevich in 
1960 and is available in an English translation 
[1960 a & b]. The theory is based on the Boltzmann 
equation for the energy distribution function of the 
electrons, thus following the classical examples of 
H. A. Lorentz and F. B. Pidduck. The authors 
point out that "to disclose the physical picture' ' the 
simpler "elementary theory' ' is "convenient and 
useful." Our present account of the "elementary 
theory" confirms this opinion. 

Similar discussions based on Boltzmann equations 
for the electronic distribution functions are con- 
tained in reports by Megill [1961], Carleton and 
Megill [1962], and Molmud et al. [1962]. These are 
directed towards special problems of considerable 
interest such as, respectively, the generation of arti- 
ficial airglows and a method to reduce the electron 
density in the P-region, by means of powerful radio 
transmitters on the ground. In all these discussions 
the use of simplifying approximations was found 
necessary. 

An interesting new way of using wave interaction 
for studying the ionosphere has been proposed inde- 



pendently by Vilenski [1962] and Ferraro, Lee, and 
Weisbrod [1963]. This depends on measurement of 
the phase changes in the wanted wave that are caused 
by the disturbing wave. 13 

A remarkable new example of the interaction of 
combination tones of two high-frequency waves with 
a third (wanted) wave, similar to that mentioned 
under prediction (11) in sections 1 and 5, has been 
reported by Cutolo [1962]. 

The waves concerned and our theoretical processes 
were as follows: 

(1) Disturbing wave: A carrier wave of frequency 
cc/2tt between 50 and 75 Mc/s was modulated at a 
frequency of Q/2ir between 1000 and 1300 kc/s which 
lies near the local gyrofrequency of 1200 kc/s; the 
peak power was 80 kw and the modulation depth 
m about 42 percent. Thus the (Yagi) aerial used 
simultaneously radiated three waves of the angular 
frequencies cu + 12, co, and co— 12. Two pairs of these 
produced in the ionosphere electric currents which 
contained components with frequencies equal to the 
combination tones (co+12)— w and co— (co— 12), i.e., 
equal to the gyrofrequency 12. These ionospheric 
current components generated secondary waves of 
the frequency 12, which in turn caused gyroresonance 
variations of the local collision frequency v. 

(2) The wanted waves were selected so as to be 
reflected at or above the 90 km level in the ionosphere 
in accordance with the requirements for gyro inter- 
action given in our section 4.9. 

Cutolo 's experiments demonstrated that in this 
way modulations of depths from 5 to 20 percent 
could be impressed on the wanted waves near the 90 
km level of the ionosphere when the modulation 
frequency approximated to the gyrofrequency. 

Lastly, we may consider the waves that could be 
generated in the interplanetary plasma by the large 
negative charge on the sun, which has been postulated 
in recent years to account simultaneously for a num- 
ber of important phenomena [Bailey, I960]. The 
evidence in support of this hypothesis has recently 
been strengthened by the fact that three predictions 
(based on this hypothesis) that were published in 
1960 have now been verified by the data on the inter- 
planetary magnetic field obtained by means of the 
space probes Pioneer 5, Explorer 10, Mariner 2, and 
Explorer 12 [Bailey, 1963]. 

The electrostatic field E s due to the charge —Q s 
on the sun is radial to the sun and has the magnitude 

E s =fQ s /r 2 e.s.u. 

at the radial distance r, where / is a number between 
and 1 which represents the partial screening of the 
solar charge by the interplanetary plasma clouds. 

Any fluctuations in the emission by the sun of these 
plasma clouds will cause corresponding fluctuations 
E in the field E s . In certain frequency bands the 
radial components of E are particularly liable to 



13 Ferraro and his coauthors mistakenly attribute the origin of the theory of 
wave interaction to its discussion in 1948 instead of to the correct date 1934. 
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amplification by the mean static field E s , since they 
may be regarded as forming longitudinal wave modes. 

These longitudinal waves may be sufficiently 
strong to produce, or help produce, certain notable 
nonlinear effects in the earth's ionosphere, which 
are at present attributed to other causes. For 
example, some geomagnetic disturbances 14 may be 
caused by such longitudinal waves [Heirtzler, 1962]. 

This possibility deserves further study, both by 
means of instrumented satellites and theoretically 
[Bailey, 1948, 1951]. 

9. Conclusion 

The foregoing account of the theory of wave inter- 
action and other nonlinear processes in the ionosphere 
in the form developed at the University of Sydney 
shows how the ideas and methods of the Townsend 
School on the motions of electrons in gases and of the 
Appleton School on magneto-ionic wave propagation 
have led to a relatively simple mathematical theory 
which has proved fertile in predictions of new 
phenomena that have subsequently been verified by 
experiment. 

It may therefore be used, and has been used, with 
confidence as the basis of new experimental methods 
for investigating processes in the ionosphere and in 
electrical gas discharges. It has also suggested 
methods of controlling parts of the ionosphere for 
particular purposes by means of beams of radio 
waves. 

The present trend towards increasingly large 
power densities in radio beams has made this theory, 
and its mathematical elaborations, of increasing 
practical importance. 

These elaborations, which are based in part on 
Boltzmann's equation for the distribution of electron 
energies, are distinctly onerous, but they serve both 
to check the simple mathematical theory of the 
Townsend School and to study ab initio those phe- 
nomena which, like the artificial airglow, depend 
strongly on the actual distribution of the electron 
energies. On account of continual and nonuniform 
fluctuations of the ionosphere and the lack of precision 
in the necessary auxiliary laboratory data, it is 
probable that for some time to come the simple 
theory 15 will serve for most purposes quite as well 
as the elaborations. 

The present account also suggests that certain types 
of low-frequency longitudinal waves of extraterres- 
trial origin may produce some modifications of the 
ionosphere that are at present attributed to other 
causes. 
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